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Abstract 

In this article, we prove the local well-posedness, for arbitrary initial data with certain 
regularity assumptions, of the equations of a Viscoelastic Fluid of Johnson-Segalman type 
with a free surface. More general constitutive laws can be easily managed in the same way. 
The geometry is defined by a solid fixed bottom and an upper free boundary submitted 
to surface tension. The proof relies on a Lagrangian formulation. First we solve two 
intermediate problems through a fixed point using mainly ^ for the Navier-Stokes part. 
Then we solve the whole Lagrangian problem on [0, To] for Tq small enough through a 
contraction mapping. Since the Lagrangian solution is smooth, we can come back to an 
Eulerian one. 



1 Introduction 

This article deals with the equations modeling the flow of a viscoelastic fluid with a free surface. 
The fluid is assumed to be incompressible, viscous and its stress tensor contains both a viscous 
and a viscoelastic part. For the proof, the latter obeys a Johnson-Segalman constitutive law 
and more general models are studied in an appendix. The geometry is 2D, horizontally infinite 
and vertically bounded by a rigid bottom and a free boundary. 
Gravity and surface tension are the only external forces. 

Unless specified, all the articles mentionned hereafter consider viscous fluids obeying the Navier- 
Stokes equations, but not viscoelastic fluids obeying more complex laws. 

A similar problem in a bounded geometry (drop of a fluid) was dealt with in numerous articles 
by Solonnikov. He wrote an article [28] in 1977 in which he proved the local in time unique 
solvability in the Holder space c 2+a,l+a / 2 (1/2 < a < 1) (no surface tension). He also proved 
the global existence with no source term and sufficiently small initial data in [27] in the space 
W 2 ' 1 with p > n (no surface tension). More recently Shibata and Shimizu [23] improved this 
result in (= 17(0, T; W 2,q (Vt)) (~) W 1,p (0, T; L 9 (fi))) and still with no surface tension. The 
latter article also contains an interesting review and related problems. Surface tension was 
included in other articles by Solonnikov where he proved local existence in time and uniqueness 
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m 

w 2 +a ,l+a/2 with 1/2 

< a < 1 for any initial data in J2B]) and global existence for initial 
data sufficiently close to equilibrium in [27]. 

An other direction of research is the flow of a fluid down an inclined plane on which Teramoto 
proved the local in time unique solvability in 3D without surface tension in [33J and with surface 
tension in [31]. Nishida, Teramoto and Win [T9] proved the global in time unique existence in 
a periodic 2D domain and with sufficiently small initial data. In [8], Bresch and Noble derived 
a shallow water model and obtained estimates e — > (still for a periodic in x flow). 

Here, we investigate an other direction in which we must quote the pioneer article of Beale in 
1981 who proved existence and uniqueness in small time in 3D without surface tension in [6], 
using Sobolev-Slobodetskii spaces. He also proved in [7] the global existence for sufficiently small 
initial data, thanks to gravity and surface tension. Fujita-Yashima proved in 1985 the existence 
of a stationnary and a time-periodic solution in the same geometry with surface tension in PU] 
by perturbation methods. Allain adapted the articles of Beale to the 2D geometry with surface 
tension. She proved well-posedness in [1] and [3J. In [29], Sylvester had large time existence 
even without surface tension. Tani extended these results to the 3D case in [31J. In a joint 
article with Tanaka, he gave a proof of large time existence for sufficiently small initial data 
whether there is surface tension or not in [32J. In [1J, Abels generalized these results to the L q 
spaces (N < q < oo) in 2 or 3 dimensions (the velocity is in W^ 1 ). In J3tJ|, Tanaka and Tani 
proved in 2003 the local existence for arbitrary initial data and global existence for sufficiently 
small initial data of a compressible Navier-Stokes flow with heat and surface tension taken into 
account. 

Recent articles were published on the numerical resolution of the latter problem especially in 
view of the simulation of surface waves. We may quote Gutierrez and Bermejo [TS] , Audusse 
et al. [5], Guidorzi and Padula [13] and in 2009 Fang et al. [9]. 

Our main result is the well-posedness of the equations of a viscoelastic fluid with a free boundary 
for arbitrary large initial data sufficiently regular. Most of the results of the present article are 
announced in [T7J and extend those of the author's PhD thesis [T6] . 

Up to some minor modifications, more complex constitutive laws can be dealt with. Basically, 
it relies on the fact that the constitutive law improves the regularity of its source term which is 
the velocity gradient. The nonlinear terms are easy to handle because of an algebra property. 
So, provided one may have estimates of the extra stress in an algebra, the same theorem applies 
to more general viscoelastic fluids. 

Below, we give the Eulerian (Subsection 2.1), and then Lagrangian equations (Subsection 2.2). 
The Equations given, we set the spaces and operators to be inverted and give the sketch of 
the proof (Subsection 2.3). Section 3 is devoted to solving an auxiliary problem which proof is 
sketched in Subsection 3.1. Then, one solves a second nonlinear auxiliary problem specific to 
the viscoelastic fluids in Section 4. Section 5 is devoted to estimates of the error terms and we 
use all the preceding results in a fixed point in Section 6. 
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2 The equations and the sketch of the proof 



Starting from the dimensionless equations in Eulerian coordinates (Subsection 2.1), we derive 
the dimensionless equations in Lagrangian coordinates (Subsection 2.2), state the operators, 
spaces and give the sketch of the proof (Subsection 2.3). 

2.1 Eulerian equations 

The dimensioned variables and fields are tilded. After exhibiting the dimensioned system of 
equations, we make it dimensionless. 

The domain of the flow is denoted Q(t) C M 2 . Its bottom Sb is given independent of time, 
and represented by a depth function b(xi) such that at the bottom z = —b. Initially fl(i = 0) 
is denoted fl The upper boundary is a free surface Spit) at time t. At t — 0, it is denoted 
Sp and represented by a height function (. We assume that these surfaces do not cross (the 
bottom does not dry even at infinity) . A typical domain can be seen on Figure [TJ 




Figure 1: Domain of the flow. 

The functions b and ( are the initial JJ 5 / 2 + r (0 < r < 1/2) height functions of the bottom Sb 
and of the free surface Sp respectively. We assume that ( tends to zero at ±oo and b to some 
limit at ±oo. So the domain is unbounded but of finite depth in the vertical direction. If we 
denote x = (xijXz) a current point in Q (so at t = 0), then Q = {x/b(xi) < xi < ((xi)}. 
Hereafter, vectors and tensors are written in bold letters. Their components are in non-bold 
letters and with corresponding indices. We use the summation convention and indices after 
a comma designate a differentiation with respect to the variable: dx Ui = Uij- We denote v 
the velocity, p the pressure, and f the extra stress tensor due to the viscoelasticity in the full 
system: 



p (d t ~v + v. VvJ - p so iAv + Vp 
div v 

Vt 

f.n — pn + 2// so zD[v].n — aHn 
v 

v(x,i= 0) 
[ r{x,t = 0) 



div f - pgoS in Q(t) x (0, T) 
in Q{t) x (0,T), 

in Q(t) x (0,T), 
(0,T), 





2/j, P oiD[v] 

Patm n 



u (x) 

T (x) 



on S F (t) x 
on S B , 
in Q, 
in Q. 



(1) 
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In this system, p is the density of the fluid, p so i the solvent viscosity, go is the acceleration of 
gravity, A the relaxation time, p po i the polymeric viscosity, a the surface tension coefficient, 
P a tm the atmospheric pressure, D[v] the symmetric part of the velocity gradient (rate of strain 
tensor) and 

T = d~ t T + V.Vf - ffa (Vv, f) (2) 

where j B (Vv, r) = ^ (Vv T r + r Vv) + s±l ( T Vv T + Vv r) , 

is the interpolated (a e [—1, 1]) Johnson-Segalman time derivative of tensors designed to let the 
tensors remain frame- invariant. The viscoelastic fluid is supposed to have its total stress tensor 
equal to the sum of a diagonal pressure matrix, a viscous term and an extra stress tensor f. 
The constitutive equation ([I]) 3 satisfied by r describes the behavior of complex fluids that have 
a memory (see [E]). Here we write the Oldroyd B model but other models could be treated 
the same way. 

To get dimensionless equations, we define a characteristic length L (uniform in every dimension) 
and a characteristic velocity Uq. They enable to define new dimensionless variables untilded: 

x = (xi,x 2 ) = Lx = L(x 1 ,x 2 ), t=^-t, (3) 
new dimensionless fields untilded: 



v(z , i) = U v(x, t) , p(x, t) = P atm - pg Lx 2 + (p so i + H P oi) -fip(x, t) , 



(4) 



and some dimensionless numbers: 



pLU XU p po i a pg L 
Re = , ' We = "T~> £ = 1 ' a = 771 i V' So = 7 : yrr- (5) 

Hsol + l^pol J-> P^sol + Upol U {p so i + Ppol) (Psoi + p po l)U 

So we denote Re the Reynold's number, We the Weissenberg number, and a the dimensionless 
surface tension. All these changes enable to make (pQ) dimensionless: 

Re(<9 t v + v.Vv)-(l-e)Av + Vp-divT =0 infi(t)x(0,T) 
divv =0 in x (0, T) 

r + We^[v]r - 2 £ D[v] = in (fl(t) x (0, T)), 

-pn + 2(1 - e)D[v] • n + r • n- aHn + g x 2 n =0 on Sp(t) x (0, T), 

v(x, t) = on S'b, 

v(x, 0) = u (x) in Q, 

t(x, 0) = cr (x) in Q. 

This system of partial differential equations is supposed to describe, in the Eulerian coordinates, 
the flow of a viscoelastic fluid submited to surface tension and gravity. We have explicitely 
used that g a is bilinear (see (j2J)), but more complex constitutive laws without any higher order 
derivative of u and of r can be included in the proof with minor changes. This will be depicted 
in Appendix A. 



(6) 
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2.2 Lagrangian equations 

We need to define the function 

fj{.,t)-. n - n(t) 

X » rj(X,t), {l) 

which gives the location at time t of the point that used to be at X G Q at time t = 0. In small 
times, rj will be close to identity and we write 

rj(X,t)=X + r)(X,t). (8) 

In a sense that will be made clearer later, the displacement r\ is "small" . Let us then define the 
fields in the Lagrangian coordinates: 

u(X,t)=v(rj(X,t),t); q(X,t) = p(r}{X,t),t); <r(X,t) = r(rj(X,t),t), (9) 

where these fields are respectively the velocity, the pressure and the extra stress tensor (mod- 
eling the polymer). We need also to define some geometrical quantities: 

(d»j)y = d Xj rg i {X,t) =fj id ( X,ty i (I) = {dff)-\X,t); 
N(X,t) = (-C(X 1 ),l)/^l + C^;Ar=(N 1 -d T r ]2 ,N 2 + d T r ]l ), 1 ] 

where <9t = (1 + C' 2 ) ^ l ^ 2 dx 1 is the tangential derivative on Sf (applied to functions that depend 
only on Xi), N is the unit normal to Sf pointing upward, and M is a vector that will appear 
later. 

In we must transform the terms to the Lagrangian coordinates on Sf- The most difficult 
term is aH(x, t)n(x, t) which is defined on x G Spit)- We will make explicit this transformation 
and let the reader check the other terms. 

Simple differential geometry gives us the unit vector N normal to Sf- In a similar vein, the 
coordinates of a current point on Sp{t) are (X\ + r)±(Xi, C(Xi), t), C(Xi) + V2(Xi, C(i), t)). If 
we denote 

Vi,X! = dxAviiXuCiXx)^)), nn 

V2, Xl = d Xl ( V2 (x 1 x(x 1 ),t)), 1 ] 

a unit tangent vector to Spit) is 

T = (1 + V i,x 1 X'(X 1 ) + r ?2iXl )/ v / (l + Vl , Xl )* + (C'(X X ) + 
^From the Frenet-Serret formula, one may write: 

Hn i dT = d T 

where dr(.) is the tangential derivative along Sf- Since it is applied on functions only of X% 



we have <9t(-) = dx 1 (-)/ y 1 + C' 2 an d this derivative is more convenient since it is geometrical. 
We may simplify further this formula by defining: 

t(x 1 ,i) = C( fj + ^'; f) -c'(x 1 ) (12) 

1 + Vi,Xi(X 1 ,t) 
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as is done in Jl] (this article refers to [3j where the reader will find some further details). This 
definition enables to write 



V(l +Vi,x 1 ) 2 + (C'(^i) + V2,x 1 ) 2 U $ ( A i.*) + C / 

(13) 

In the Lagrangian coordinates, we will need also an evolution equation for $(Jfi,t). Since 

^r7(X,t) = u(X,t), 
one finds easily the time derivative of $: 



Easy computations enable to complete the derivation of the Lagrangian system of equations 
that may also be found in classical books ([IS], [22], ...): 

' Re m t t - (1 - e)l,kMijUi,i),k + IkAk - Oij,kl k j = in ft x (0, T), 

=0 inftx(0,T), 

Oij + We ^tcrij - gL ^(^ li u ky ia kj + a ik u k ^ 

-£(ui,k£ k j + Uj,k£ki) =0 inftx(0,T), 

crijKj ~ qKi + (1 - e)(£kj u i,k + Ui u i,k)Kj+ 

+g (ax 1 )+7 l2 (X 1 ,t))K l - (14) 

a (d r ({I + ($ + C') 2 )^ 1 ( * ! w HI =0 on S F x (0,T), 



$t _(druM/ =Q on ^ x(0 ,T), 

$(t = 0) = on S F , 

u(X,0) =u (X) in ft, 

er(X,0) =cr (X) in ft, 

u =0 onS B x(0,T). 

Here, e G [0, 1[ is the dimensionless polymeric viscosity, go the acceleration of gravity, N the 
outward unit normal to Sp, and J\f = (Ni — 9 T ^ 2 , N 2 + 9 T ^i) a (non-unit) vector normal to 
5V(i) at the point Tj(Xi,£(Xi),t) convected back on ft, defined on (Xi,t) and geometrically 
defined with <9t instead of dx 1 - 

In system ( Tlijl . there are equations inside ft ( f|T4l) i to ffLljU ). equations at the free boundary 
of ft defined on X 1 G R ((p3D 4 and flUDs), initial conditions ( ([13)1 w to <|I3jl n ) and a Dirichlet 
condition on the bottom fjT4|) o. 

/,From now on, the physical domain is the initial one and we consider only the Lagrangian 
formulation. 
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2.3 Operators, spaces and sketch of the proof 



2.3.1 The operators 

Let us remember that for short time, rj(X,t) = X + r](X,t) ~ X. So the displacement rj will 
be small (in small time) in a sense to be defined. For the same reason, if we define £: 

((dry)- 1 =) I = Id + £, 

we see that £ = (Id + d^) -1 — Id will be small. 

If we define the operator P(£, u, q, 0, cr) as the left-hand side of (THj) for equations <HM i to (ITljk 
and the initial conditions (TTilW to ffT4l) a. then ffT4j) amounts to solving 



P(£,u,g, 



(0,0,0,0,0,u , o-o) 



for u vanishing on S B (see ([HJlg) and $(t = 0) = (see ffT4L). 

The orders of magnitude of various terms must now be identified so as to see (Cbil) as a pertur- 
bation of an invertible system. 

In system fTH|) . there are some source terms : gravity and the initial curvature that are not 
small even for small times. They are of zeroth order and will be put in P(0, 0, 0, 0, 0). 
All the terms containing at least one £ are stored in E (for "Error"). The occurrence of a £ 
will enforce smallness. 

Proceeding as in [1], we collect the remaining terms in a new operator acting on (u, q, 0, cr) with 
the zeroth order in £ (we put £ in E). Moreover, we keep only the linear in $ operator from 
(TH|) and put the h.o.t. in $ terms in E 5 (since $ is initially vanishing these nonlinear terms will 
be small in small time). We will denote Pi(u, q, 0, cr) this new operator with a minor change 
= (1 + C' 2 )" 1 ^- So an auxiliary problem to be solved is Pi(u, q, 0, cr) = (f, a, m, g, k, u , cr ) 
and writes: 



Red t u - (1 - e)Au + Vq - 


diver 


= f 


in Q 


X 


(0,T), 


div u 




= a 


in Q 


X 


(0,T), 


cr + We (d t cr - # a (Vu, cr)) 


-2eD[u] 


= m 


in Q 


X 


(o,n 


<r-N-gN + 2(l -e)D[u] 


•N-a<9 T (0N) 


= 9 


on Sf 


X 


(o,n 


<jH - (<9 T u) • N 




= k 


on Sf 


X 


(o,n 


u(t = 0) 




= MX) 






in f2, 


cr{t = 0) 




= cr (X) 






in f2, 



(15) 



for initially vanishing 0, and u = on Sb- It is a viscoelastic Stokes flow. Notice that passing 
from (fl4|) to (fT5j) . the gravity term (yf (C(^i) + ^2) A/" and the zeroth order surface tension 
disappear because they are in the zeroth order term P(0, 0, 0, 0, 0). In the Navier-Stokes case 
Pi is linear (Stokes) but because of the nonlinearity of g a , our operator Pi is nonlinear. 
With the preceding definitions, (CHI) is equivalent to 



P(£, u, q, 0, cr) 



P(0,0,0,0,0) + P 1 (u,g, 
(0,0,0,0,0,u ,<t ), 



+ P(£,u,g, 0,<r) 



(16) 
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for a velocity vanishing on the bottom and $(t = 0) = 0(t = 0) = 0. In the remaining, we 
will need also to define the nonlinear auxiliary operator P2 which gives the nonlinear behaviour 
close to (ui, gi, 0i, <Ti): 



P 2 [ui,<ri](u,g,0, 



Pi(ui + u, q 1 + q, 0! + 0, a 1 + cr) - Pi(u 1; q 1 , 0i, cr x ) 
( Re d t u - (1 - e)Au+ Vq - diver \ 
divu 

cr + We (<9 t cr - # a (Vu, cr) - # a (Vui, cr)- 

9«(Vu, <ti)) -2eD[u] 
a ■ N - qN + 2(1 - e)D[u] • N - a<9 T (0N) 
0t — <9tu ■ N 
u(t = 0) 

V = o) / 



(17) 



for u vanishing on 5b and 0(t = 0) = 0. In the Navier-Stokes case P 2 = Pi because Pi is linear. 



2.3.2 The functionnal spaces 

Following J.T. Beale [B] and G. Allain [4j, we define anisotropic Sobolev-Slobodetskii spaces for 
any s: 

K S {Q x (0, T)) = L 2 (0, T; H S {Q)) f| PP /2 (0, T; L 2 (ft)), 
^(SV x (0, T)) = L 2 (0, T; H°(S F )) f| PP /2 (0, T; L 2 (5 F )), 

where PP is the classical Sobolev space whose norm is denoted | . \ s . We denote | . \k« the 
norm of the space K s . The properties of these spaces, namely embedding and trace theorems, 
can be found for instance in [IB] . 

One may define the space of (u, q, 0, cr) for < r < 1/2: 

X r T (nx (0,T)) = { (u,g,0,<r)/ 

u G K r+2 and u = 0on5 B x (0, T); 

Vq G K r and q\s F G H^+l (0, T; L 2 (S F )); (18) 
<9 T G L 2 (0,T;PP+§(,S F )) ; 0, G iT+l^ x (O,T));0(O) = 0; 
(7 6^(0,7;^+^)) }. 

This space XJj, is basically a regularity space with the condition at the bottom for u and the 
initial condition of 0. The (0, T; P 1+r (fi)) regularity on cr is needed to have an algebra 
so as to manage the nonlinear terms of a wide variety of constitutive equations. 
Indeed, the Vucr terms of the constitutive law must be estimated in L 2 (0, T; H 1+r ) f] P^~(0, T; L 2 )\ 
For L 2 (0, T; P 1+r ), if we use the L 2 (0, T; P 1+r ) for Vu, we need a L°°(0, T; P 1+r ) estimate for 
cr. For P^(0, T; Z 2 ), if we use the P^(0, T; L 2 ) estimate for Vu, we need a P^(0, T; L°°) 
for cr. So estimates of cr in ^^(OjT; H 1+r ) are natural. 
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The image space Yf is 



Yf(Qx(0,T)) = { (f,a,m,g,k,uo,tr )/ 
f G K r {Q x (0,T)), 

a G L\0,T;H^\Q))C]H^(0,T, oH' 1 ^)), 
meK r+1 {Qx(0,T)), liyj 
g,keK r+1 HS F x(0,T)), 
u ,cr GiP +1 (Q)}, 

where oP _1 (^) is the dual space of °H 1 = {p G H 1 jp = sur Sp}. We will need to keep the 
space #§ +1 (0,T; qH- 1 ^)) so as to have a lift field u G H% +1 (0, T, L 2 (f2)). 



2.3.3 Main result and sketch of the proof 

Our main result of well-posedness is the following: 

Theorem 2.1. Let < r < 1/2, and the initial height functions of the free boundary ( and of 
the bottom b be such that (—b — limb) G H r+5 ^ 2 (IR) . Let also (uo,cr ) be in H r+1 x Hl^l and 

±oo y 

the compatibility conditions 

div Uq = in Q, 
u = on Sb, 
<t .N.T + (1 - e)(D[uo]).N.T - ad T (T) .N.T = 0, 

where N is the unit normal and T is the unit tangent vector, be satisfied. Then there exists 
To > depending on the data ; r, Q, uq, cr , (, b, We, e, a and there exists a unique (u, q, 0, a) G 
solution to the Lagrangian system [Lfy . Under the same hypothesis, the Eulerian system 
(TIP admits a solution with r\ G H 1 ^, T ; H 2+r (tt)) C] F 2+ i(0, T ; L 2 (fi)). The solution depends 
continuously on the initial conditions provided they are in a bounded subset of (H r+1 (Q)) 2 . 

Remark 2.2. One could even prove that the solution depends continuously on the parameters 
(e, Re, We, a) in the domain [0, 1) xl +2 x [—1, 1]. To that end, we should track the occurrence of 
the parameters in the constants. Since the constants in the estimates depend only polynomialy 
on the parameters we can have the same estimates on any parameter in an open neighborhood 
of the given parameter. This continuity is very unlikely to be extendable to the whole domain 
at least since e = 1 would not have the regularizing effect of the Navier-Stokes equation. 

In a first part (Section 3) we solve (|T5|) also written Pi(u, q, (ft, cr) = (f , a, m, g, k, Uq, <To). 
For that purpose, we first introduce a reduced auxilliary problem P 2 [ui, <Ti](u, 9, 0, cr) = 
(f , 0, m, 0, 0, 0, 0) with zero initial conditions (lifted in the Ui,(Xi) that we solve through a 
fixed point between the Navier-Stokes part (solved in and the constitutive equation. This 
step needs uniform (Subsection 13.31) and contraction (Subsection I3.4|) estimates. Using this 
auxilliary problem, we solve ffT5|) and so invert Pi in Subsection 13.61 In Subsection 13 . 71 we show 
the continuous dependence with respect to the right-hand side (rhs) of fTl5l) which comprises 
the initial conditions. To show that the final mapping P{~ is Lipschitz, we have to check that 
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this inverted operator is bounded (independently of T < T ). 



In a second part (Section 4), for given (ui,0i) inabounded ball of K 2+r xK r+1 f] L°°(0, T ; H 1+r ), 
we prove there exists a unique solution in X£ (= Xj, with zero initial conditions) to the prob- 
lem P 2 [ui,(Ti](u,q,(j),(T) = (f , a, m, g, k, 0, 0) and prove boundedness of P 2 [u 1 ,cr 1 ]~ 1 with a 
constant independent of (ui, cr\) in a given ball. 

In the third part of the proof (Section 5), we derive some estimates on the "error" terms E 
insuring that they are small and contractant for T small enough. 

Last, in Section 6, we gather these results to build a functional F involving Pi, P 2 and E that 
is contracting. The regularity found ensures then the existence and uniqueness of a solution to 
the eulerian problem fTT4"j) . This completes the proof. 

3 Inversion of the operator Pi 

In the present section, we prove the following Theorem of well-posedness: 
Theorem 3.1. Let < r < 1/2, B > be given and the compatibility conditions: 



be satisfied. There exists a time To such that ifT< Tq, for any (f, a, m, g, k, Uq, <t ) in the ball 
in Yj, of radius B denoted B Y f{0,B), there exists a unique (u,p, 0, cr) G X£ solution of (TI3j). 
Moreover the (nonlinear) operator (u,p, (ft, cr) = Pf (f, a, m, g, k, u , cr ) is such that 



where C = C(e, We, B, T , Re) is a constant that does not depend onT < T Q . 

Proof. To prove this Theorem, we need to solve a restricted problem P2[ui, <Ti](u, q, <fi, cr) = 
(f , 0, m, 0, 0, 0, 0) with vanishing initial conditions in which ui and <j\ lift the initial condi- 
tions. As a consequence, notice that the continuity with respect to the initial conditions means 
continuity also with respect to those fields. Subsections 3.1 to 3.5 are devoted to this, while Sub- 
section [223 proves existence of Pf 1 and Subsection 13.71 its boundedness (j20|) . Local continuity 
with respect to the initial conditions is a mere consequence of (120|) . 

3.1 A reduced second auxiliary problem 



Let (ui,o-i) be in K r+2 x K r+1 f| L°°(0, T ; H 1+r ). To prove Theorem EU we need to solve 



div Uo = in VL, 
u = on Sb, 
(o- 0i)i ).N.T + (1 - e)(D[u ]).N.T - ad T (T) .N.T = 0, 



| u-u', p-p',<f)-<f)', cr -cr' \ x ^< 
C I f — f, a — a', m — m', g — g',k — k' 



, U — Uq, (Tq — (Tq |yr, 



(20) 



P 2 [u 1 ,cr 1 )(u,q,^,(T) = (f,0,m,0,0,0,0). 



(21) 
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for (u, q, $, a) G X r T . 

The resolution of the non-reduced problem P 2 [ui, <Xi](u, g, 0, er) = (f , a, m, g, k, 0, 0) is post- 
poned to Section H] and will use the resolution of (J2"T|) . 

To prove Theorem 13.11 in a first step, we solve (121"]) in Subsections 13.31 and 13.41 In Subsection 
I3.6[ we invert Pi and in Subsection 13.71 we prove that its inverse is bounded (12"U|) . 
Let us begin by a very brief sketch of the proof, then a more precise justification of the required 
estimates. We report the beginning of the proof to Subsection 13.31 after recalling some lemmas 
in Subsection 13.21 

In the system (l2~Ti) . where P 2 is defined in ( fT7|) . we will split the equations with the Stokes/Navier- 
Stokes part on the one hand and the constitutive law on the other hand. This splitting is done 
by numerous authors (see e.g. Guillope-Saut [12] and Renardy J23])- Then, given u n ,q n ,(p n , 
we will find an extra stress er n+1 . This extra stress may be carried forward in the Stokes-like 
system solved by [I] and will provide u n+1 , q n+1 , <fi n+1 . 

Then we will prove that this sequence u", g n , <fr n , & n is a Cauchy sequence in a Banach space. 

Let us be more precise below on the required estimates for the full proof. Let u° = 0, q° = 
0, 0° = 0,(7° = 0. We assume we know (u™, q n , n , cr") G X^ o and we look for er n+1 such that: 

J o-n+i + We (d t v n+1 - <7„(Vu n , <r™ +1 ) - ^(Vu 1; cr n+1 ) - g a (Vu n , trj) = 2sT>[u n ] + m 
\ cr n+1 (t = 0) = 0. 

(22) 

Such an equation is a simple ODE in the space K 1+r with no loss of regularity thanks to the 
Lagrangian simplification. Even, there is a gain of regularity of cr n+1 with respect to Vu" 
because there is no time derivative on Vu" and no gradient on a n+1 . The nonlinearity will 
be managed with uniform bounds on I <r n I i+r , . Thus, the velocity will occur only to 

& 1 1 H^2T (0,T;H 1 + r ) ' J J 

the source term Vu n . This enables to gain one level of time regularity and prove estimates 
of er n+1 in H 1 ^ 1 (0, T; H 1+r ). Indeed, since Vu b a x , Vu n , m are in L 2 (0, T ; H 1+r ), we have 
(T n+1 E H l (<d,T -H 1+r ) C H^^.T^H 1 ^). Notice that this is true until time T . This 
gain will enable to prove that (u ra , q n ,(p n , cr n ) 1— > (u n+1 , q n+1 , ra+1 , cr n+1 ) is a contraction for T 
sufficiently small in X^. 

Then we look for (u n+1 , q n+1 , <p n+1 , 0) G X£ solution of the Stokes like problem: 

f + div<r n+1 in tt x (0,T), 
in x (0,T), 
-<r n+1 • N on S F x (0,T), 
on S F x (0,T), (23) 
on Sp, 
in tt, 
on 5 B Vt. 

Theorem 4.1 of [4] in the case f = f + divcr n+1 ,a = 0, g — —er n+1 ■ N, k = 0, Uo = states 
that for given rhs and To (even though it is not obvious in the way it is written, one may state 



Re d t u n+1 - (1 - e)Au" +1 + Vg n+1 
divu n+1 = 
-<f +1 N + 2(1 - e)D[u n+1 ] ■ N - ct«9 T (0 n+1 N) = 
< - (d T u n+1 ) ■ N 

n+1 (O) 
u" +1 (0) 

u n+l 

V 
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also this result thanks to the regularizing effect of the Stokes operator), there exists a unique 
(u" +1 ,g n+1 ,f +1 ) in X r To if f + div<r n+1 G K r (tt x (0,T )) and a n+1 ■ N G K r+ l(S F x (0,T )). 
Thanks to Theorem 4.1 of [4], we have 



u n+1 , q n+ \ n+1 , \xr<C\f + divcr n+1 +C \ a n+1 



K r+ i(S F x(0,T)) , 



where C will denote a constant independent of T < T until the end of this article unless 
otherwise mentioned. From the assumptions and the definition of cr n+1 , we have that f + 
divcr™ +1 G K r . Still thanks to the fact that <r n+1 G H\0, T; H 1+r ) we know that cr n+1 ■ N G 
^ +1 / 2 (5 F x(0,T )). So 

I ,,n+l „n+l j.n+1 n I << H \ f \ J~C* I sr n + l I 

I u ,g ,0 ,U | Xf < C \t \ K r +G | cr l i2(0iT . H1+r)nif 5(o,T;ifi)nifi + ^(o,T ; ^) • 

So, for (u n+1 , g ra+1 , n+1 , 0), we will need uniform estimates of <x™ +1 in various spaces. 

For (0, 0, 0, er n+1 ), we will even need to estimate cr n+l in H^~(Q, T; H 1+r ) which is even 

stronger. So we will prove only the latter. 

Concerning convergence, since the Stokes-like system is linear (it is not the case of our consti- 
tutive equation !), we even have thanks to the same Theorem 4.1 of [3]: 

I u n+1 - u n , q n+l - q n , <f +l - 4> n , 0\ X r<C\ div(cr" +1 - a 71 ) \ K r +C I <x™ +1 - a n I r+ 1 

Thanks to the continuity, injection and trace theorems found in [TSj (Chapter 4) or in [6] 
(Lemma 2.1): 



div(<T n + 1 — <T n ) |L2( 0iT;ff r)< C I <7 n+1 — CT™ |L2( 0jT;/i -l + r) 

div(<7 -CT ) | h5(0 ,T;L 2 )- 6 I _ 



r2(0,T;L 2 )— I l/f^T;/*" 1 )' 
_n+l I ^ C* I <X n ~^ 0" n I 



So 



I u n+1 - u n , q n+1 - q n , n+1 - n , | X r < C I <r n+1 - <r n I 1X „ r , lx „ -.i i . (24) 

To prove that the (u n+1 , q n+1 , <f) n+1 , 0) part of our sequence is of Cauchy type in a Banach 
space (and we will always need C independent of T < To), it suffices to prove thanks to the 
interpolation Lemma [3.51 

I <T n+1 - I 



L 2 (0,T;H 1 +'-) 

l + r 

m~(0,T;L 2 ) 
H%(0,T;H l ) 
H2+i(0,T;H2) 

where e' is a positive constant not depending on To nor 
In addition, for (0, 0, 0, cr n+1 ), we need to prove also 



_ra+l _n i 

(T — (T l + r 



CT n+1 



<J — (T r , 1 



< CT e ' 


| u" 


- u n 


-1 




< CT e ' 


u n 


- u n - 


-1 




< CT e ' 


| u n 


- u n 


-1 




< CT e ' 


| u n 


- u n 


-1 


\k<~+ 2 



n. 
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(J — <T \ l + r L/i U — U 

1 l ffi"(0,T;I 1 +'')- 1 

which contains all the others. So we will prove only the latter. 



K r + 2 , 



As a first step for all the preceding estimates, we need uniform in n and in T < To estimates 
so as to manage nonlinearities: 

3T,V,S/\ (T n \ m+ r (t) <S Vn, < t < T (25) 

\u n ,q n ,<f) n ,0\xz <V Vn. (26) 

This will be done in Subsection 13.31 after some lemmas stated in Subsection 13.21 

As a second step, we prove in Subsection 13.41 that the sequence is of Cauchy type thanks to the 

fact that C is independent of T < To and to the T e term: 

• | cr"' +] - cr n \ H ^z {0>T . H1+r) < CT e | <x n+1 - cr n \m(o,T-,m+r), (27) 

Vu n_1 \L 2 (0,T;H 1 + r ) (28) 

l_1 \ K r+*, (29) 

Among these inequalities, (1291) is obvious. The inequality (|27|) will be a simple consequence of 
Lemma [3.41 proved hereafter that relies on (cr n+1 — cr n )(t = 0) = 0. The last inequality will be 
proved in Subsection 13.41 After lifting some rhs, this will solve (12"T|) . 



< 


CT € ' 


a n+l_ 


< 


CT £ ' 


Vu"- 


< 


CT € ' 


u n -u' 



Coming back to the proof of Theorem 13 .1[ we will lift the initial conditions with ui, q%, 0i, (J\ in 
Section ESI We will also justify that a 1 e K 1+r f| L°°(0, T; H 1+r ) to apply Theorem ITU Then 
the field (ui+u, qi + q, <pi+(p, er 1 -\-er) G Xj, will solve ffl5|) with non-vanishing initial conditions. 

In Subsection 13 .7\ we prove uniqueness of the solution of (IT5|) . The estimates enable to prove 
that the inverse of Pi is bounded independently of T < To. This will complete the proof of 
Theorem 13.11 



3.2 Some useful lemmas 

Lemma 3.2. Let X be a Hilbert space, < s < 2, such that s — \ is not integer. 

There exists a bounded extension operator from ju G H s (0, T; X), c^u(0) = 0, 0<A;<s — ^| 

in H S (M. + ; X) . The constant C of its boundedness does not depend on T < Tq. 

The proof of this fundamental lemma based on reflections is classical and can be found either 
in [T8J p. 13 or in [6 J (Lemma 2.2). 

Remark 3.3. The initial vanishing conditions (even up to some derivatives) are needed what- 
ever the extension operator (see remark \57fy . This is the reason why we have to solve the system 
with vanishing initial conditions. For such functions the constants do not depend on T < Tq. 
Then we lift the initial conditions to solve the full system. 
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Lemma 3.4. Let X be a Hilbert space, 1/2 < r < 1, < s < r and s / ^. Iff G iF(0,T;X) 
and /(O) = ; i/ien ; 

I / \h s (0,T;X)< CT T S I / \H r (Q,T;X) ■ 

The constant C does not depend on T <T Q . 

This lemma is proved in [3] (Lemma 5.3). The proof is based on extension and interpolation. 
We need to have /(0) = to have C independent of T < T . Indeed, \e~ t \L2(o,T)/\ e ~ t \H 1 (o,T) 
can be bounded by a constant that does even not depend on T < T . 

We will repeatedly need the following lemma about identity operator from K r to the space 
H p (0,T; H r - 2p (tt)): 

Lemma 3.5. Suppose < r < 4 and p < £. 

(i) The identity extends to a bounded operator K r (Vt x (0,T)) -> fP(0,T; H r - 2p (9)). 

(ii) If r is not an odd integer, the restriction of this operator to the subspace: 

e K r (Q x (0, T)) / d*v(0) = 0, VA; < fc < 
bounded independently ofT<T . 

Proof. 

(i) From the definition of K r , identity can be defined from K r to L 2 (0,T; H r ) and to 
H%(0, T; L 2 ). It is bounded. So identity can be defined on the interpolate space be- 
tween K r and K r (which is K r !) to: 

[L 2 (0, T; H r ), Hi (0, T; L 2 )] = H* (0, T; H r ^), 
from Proposition 4,2.3 of [TS]. It is the announced result with p = 9r/2. 

(ii) Lemma [3.21 guarantees that the extension to t > T is bounded independently of T < To. 
Then one may apply the same proof as above on [0, +oo] x Q. The norms are bounded 
independent of T < T . 

□ 



3.3 Uniform estimates 



We want to prove (1251 1261) by induction. Indeed we will prove a more general result. The 
inductive hypothesis writes for (u n , q n , (fi n , er n ) defined by (j2"2~j |23|) : 

I <r n \m + r (t) < S, 



3T , V < 1, S positive such that Wn < 



u n , q n , 







x. 



<V, 



|L 2 (0,T ;H 1 + r ) 



n< V, 



< 1, 



C I VUi |l 2 (0,T ;// 1 +'')< 1) 
Cy/%(V+ | VUi |l2(0,T ;H 1 +'-), 

C(y+ I m |L2(o jTo;if i+r)) < 5, 

C I f l^5(0,To)< 

[ CT^'(V+ | m | i2(0) T 0i ^)) < V/5. 



(30) 
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Notice that the five last inequalities do not depend on n. So it suffices to determine a To at 
n = to satisfy these properties. Notice also that the L°°(0, T; H l+r ) estimate on cr n+1 cannot 
be omited in the hypothesis because of the nonlinear terms in the constitutive equation which 
will require even more regularity. Notice also that (13"0"jU is included in (JSP]) -?, but more useful. 

Obviously the property at n = is true if To is small enough. 

Let us assume property (130]) be true for n and let us try to deduce that it is true for n + 1. We 
must be careful that the constants do depend neither on T < T nor on n. 

In a first step, we want to prove (1251) or (l3"Uj) i from the H 1+r scalar product of (I2"2"j) with cr n+1 . 
Thanks to the fact that g a is only quadratic, it provides (with a C independent of T < To and 
of n)\ 



-n+l |2 , We Cl I n+l |2 ^ ^ 



2e I Vu" \ 1+r + I m |i +r + 



+ CWe(| Vu™ | 1+r | | 1+r + | V Ul | 1+r | | 1+r + | Vu™ | 1+r | cr | | 1+r ) 
and so 

(1 - CWe(\ Vu" \ l+r + | Vm \ 1+r )) | <r" +1 \ 2 1+r +^j- t I cr n+1 !?,.< 



_n+l I 
CT |l + r; 



(31) 



ll+r- 



C [(2e + CWe | (T l \ 1+r ) \ Vu" | 1+r + | m | 1+r ] | a 



n+l 



\l+r 



If we choose not to simplify | cr n+1 | 1+r , Holder inequality enables to state: 



_d_ 
dt 



( (1/2 - CWefl Vu" | 1+r + | Vui |i +r )) 



2 

e 



We 



\ 



.n+l 1 2 



V 



< 



/ 



C( £ , < r 1 ,We)(| Vu" |? +r + | m |? +r ) e 



;i/2-CWe(| Vu" |i +r + | Vm |i +r )) 
We 



ds 



Thanks to the assumption that cr 1 E L°°(0, T; H 1+r ) and that a n+1 (t = 0) = 0, one may 
deduce: 



1 (1/2 - CWe(| Vu" | 1+r + | Vui | 1+r )) 
We 



t -2 

(e "'s 



(C | Vu" \ 2 1+r + | m \i +r ))ds. 



o 



We have to estimate the term in the exponential. Since s < t and thanks to ([3"01 for n: 



-2 



l/2-CWe(| Vu' 



l+r 



Vui 



l+r) 



We 



ds <C I | Vu" |i +r + | Vui |i +r ds' 



< Cy/T (j Vu" |l 2 (0,T;H 1 +'') + | Vui |l 2 (0,T;// 1 +'')) < 1- 



This enables to state: 



n+l 



l+r (*) < C (| Vu" |l2( 0i T;H 1 +'-) + | m \l 2 {0,T;H 1 +^)) < C (V + \ HI |L2( 0iTo;H i+r)) , (32) 
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and we deduce the announced (j25j) or (1301) i for suitable 5 not depending on n. Moreover we 
have: 

< CT e (V+ I m \ L 2^ T , m+r) ) < S. (33) 

Remark 3.6. The proof above relies on the fact that the C found at n+ 1 th iteration does not 
depend on the C assumed at the n th iteration. 



In a second step, so as to have uniform in n estimates on | u n \x 2 + r (denoted (130|) 9 and (150)1 q). 
we need more regular in time estimates on er. More precisely we need | cr n+1 \m(Q,T;m+ r )- 
We take back (T22l) . isolate the time derivative and take the H 1+r (Q) norm: 

IS^Ii+r (t) < ^ | Vu" | 1+r +^ | m | 1+r +^ | \ 1+r + 

+ C(\ Vu" | 1+r | a n+1 \ 1+r + | V Ul \ 1+r \ a n+1 \ 1+r + | Vu* 1 \ 1+r \ n | 1+r ) . 

One may then take the L\ norm and it writes thanks to (l30|) . (1321) and (1551) : 

|a t <T n+1 | i2(0T;H1+r) < C(e,We) (| Vu" | L 2(o,T;tfi+'') + I m \ L ^(o,T;m+^)) + 

+ C|cr n+1 | L 2( 0iT . H l + r) + C\ Vu n | L 2( 0iT;H l + r)|cT n+1 | L o O ( 0iT . j yl + r) + 
+ C| VUl | i 2( 0iT; ^l + r) ICT™" 1 " 1 | Lo o( 0)T .^l + r) + 
+ C | VU™ | L 2 ( 0iT .^l + r ) | O"! | L oo ( 0iT; ^l + r) . 



It is then easy to use (1331) to prove 

I c n+ \m(o,T;m+>~)< C(e, We, ui, cr 1; 5)(| Vu n lia^/rn-*-) + | m \L 2 {o,T,m+ r ))- (34) 

This estimate is even stronger than the estimates required in L 2 (0, T; H 1+r ), Hi (0, T; i? 1 ), and 
#5+4(0,T;#2). Because cr" +1 (t = 0) = we may use Lemma 13.41 and Lemma 13.51 which imply 
that the constant C below is independent of T < T : 

I ^(O.T;// 1 )- | <T n+1 Ih^O.T;// 1 )^ CT € (| Vu" \l 2 {0,T,H 1+t ) + I m I L 2 (0,T,H 1 + r ) ) • (35) 

Similarly, Lemma [3.41 and Lemma [3.51 enable to get: 

I <T n+1 I r i 1 , h< CT" 6 I (T n+1 I , , 1 < CT 6 (I VU™ \j2(n T.TTl+r\ + I III \r2/n r.ffl+rO 

(36) 

We are in position to get estimates on | u n+1 \x 2 + r independently of n and in particular 
(150)1 o. Indeed we use the Theorem 4.1 of G. Allain [4] with # = -<x™ +1 • N, k = u = a = 0, 
f = f +div <t™ +1 to find solutions to (1231) . Of course we use also the already proved estimates (1551 
[35], [36]) together with classical trace theorems and the induction hypothesis ([30]) . If < r < 1/2, 
it enables: 



,,ri+l „n+l 



m+1 



,0| x ,<C|f + div<T« +1 | K 



< C | f 
^ C f Ik 1 " 

< V. 



K r +C I 0" n+1 



L 2 (0,T;_ff 1 +'")n^ 2 (0,T;i ; / 1 ) 



,n+l 



-CT e '(V+ I m \ L *(o,T;Hi+r)) + CT e '(V+ I m |£2 (0j7W i+r))+ 

+ CT e \V+ | m |x2(o,T;fl-i+'-)) + CT e '(V+ | m |z2 (0)T;iJ -i+r(n)) 



l K r+ ?(S F x(0,T)) 

L 2 (0,T;H 1 + r (n))n-H'5 + 3(0,T;_ff^(n)) 



(37) 
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In the inductive hypothesis we assumed V < 1 and S whatsoever. Should we take a smaller To 
(the constants C above depend on e, V, S, We, Re, cr 1; T but not on T < T nor on n) we may 
derive (130]) 9 and (l30lU for n + 1 assuming the inductive hypothesis for n. 
Since the inequalities flSUL to f[3"Uj) a even do not depend on n, they are satisfied from the fact 
that they are satisfied at n — 0. So the inductive proof of (13T)j) is complete and (121)1 1251) are 
also proved. 

The case of more general constitutive laws is studied in Appendix A. 



3.4 Convergence of the sequence 

We intend to prove that the sequence (u n , q n , <p n , <r n ) is a Cauchy sequence in the Banach space 
XJp. Although we proved (I25H26P and even fl30|) . we will only use: 

I \L x (0,T;H 1+r )^~: S \/fl. 

We still need the estimate (128]) that will be proved hereafter. Notice that our domain is 

unboudned. This would make compactness argument very difficult to use. 

In order to prove convergence, we take the difference of (T2"2"j) at the n th and n + 1 th iteration : 



_.n+l _n 
<J — (7 



We ( d t {a- n+1 - <r n ) - g(Vu n , cr n+1 - cr n ) - g(Vu n - Vu n ~\ cr n ) - 
p(Vui, cr n+1 - a n ) - g(Vu n - Vu"" 1 , tr{) = 2e~D[u n - u^ 1 ]. 



We want to prove (128]) . The uniform estimates (I25H26I) and assumptions on er 1; Vu x enable to 
simplify the scalar product in H 1+r of the equation (1551) with cr n+1 — er": 



+CWe 



_n 12 j_ We ^ I „n+\ n |2 _ I T7„n 

ll+r+^-^l' 7 ll+r< 2£ | VU - 

V-,-. n I I |2 i I Y7,, ra T7n ra — 1 I I ~n I I ^"+1 I _L 

u | 1+r | cr T — er | 1+r + | VU — VU | 1+r | CT | 1+r | (T T - IT | 1+r + 



rt-1 1 

|l- 


-r 




— er 


^ |l- 


-r 


CT n+l 


— er 


r 1 ^ 


-1 


- <T n 


l+r] 



+ I V Ul \ 1+r \ (T n+1 -o- n \ 2 1+r + \ Vu n - Vu"- 1 | 1+r I a x \ x 

(1/2 - CWe(| Vu" | 1+r + | Vu x | 1+r )) | cr^ 1 - cr- |? +r + W§ A | er" +1 - cr" \ 2 1+r < 

C{e, We, V, 5, cr x ) | Vu n - Vu n ^ \\ +r . 

As in the proof of (l33|) from (l3T|) one may use (1301) to have: 

| cr" +1 - cr™ \ l+r (t) < C(e, We, V, S, cr,) | Vu" - Vu^ 1 | LW; ^) • (39) 

We want to prove (|2"%1) . Since we already know that | er n |l°°(o,T;_h" 1 +'-)< S, and | er 1 |l°°(o,t ;J3" 1 +'-) 
and I Vui I , , i+r „ are in a bounded ball (here we need this to have continuity 

for bounded initial conditions because ui, a\ will depend on them), we may use the inequality 
(l39|) and come back to (l38|) whose L 2 (0, T; H 1+r ) norm provides: 
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\d t {cT n+1 - OI^(0,T;ff^) < I " Il 2 (0,T;^+") +2^ | Vu" - Vu""' ^(o.T^+r) + 

i I T7,, n I I _n+l _n I i 

+G | VU |l 2 (0,T;H 1 +'')| CT — <T |l°°(0,T;-H" 1 +<-) + 

+C(S) | Vu" - Vu"- 1 Ir^ffi+r) + 

+ C I VU! |l 2 (0,T;_H' 1 +'')| CT n+1 — CT™ |ioo(o ) T;H 1 + r ) + 
+C | CTl |l°°(0,T;/i 1 +'')| Vu™ — Vu" 1 \L 2 (0,T;H 1 + r ), 

< CfaWeMS,^,*!) I Vu* 1 - Vu"" 1 \ LHo , T ., H i +r) . 

Notice that the constant C(e, We, V, S, u 1; o^) depends on Ui, <T\ only through a bound of their 
norm. The latter can be combined with (1391) to get the requested inequality: 



cr — cr 



H 1 (0,T;H 1 + r ) 



< C | Vu" - VU™" 1 | L 2 {0)T;H i +r) . (40) 



3.5 Passing to the limit 

Since estimates (J27j) to fl22D are proved for (0, 0, 0, cr™) and estimate (T3U) for (u n , q n , $ n , 0), the 
whole sequence (u n , q n , (fi n , er n ) is of Cauchy type, should it be for sufficiently small T . So it 
converges in X? to some (u, q, cf>, cr). 

Can we take the limit n — > +oo ? When the terms are linear it is straightforward. When they 
are of the shape Vu™ cr™ in the constitutive law, we remember that Vu™ converges in K 1+r 
and cr n in (0,T; H 1+r ). So 

| VuV - VUCT \ K l+r < C | VU™ - VU \ K l+r\ <T" \ ^ ( 1+r) + 

Notice that a priori we use the bound C of the embedding H~ (0, T) > L°°(0, T) that depends 
on T < To. But we prove further Lemma 15.91 which enables to state that the bound can be 
taken independent of T < T because (Vu™ - Vu)(t = 0) = and <r™(i = 0) = 0. So Vu"cr™ 
tends to Vucr in K 1+r . Even if more complex nonlinearities of Vu and cr had been chosen, 
they could pass to the limit provided uniform H~ (0, T; H l+r ) estimates of cr™ can be derived. 
So the most difficult step for more complex constitutive laws is the uniform in n estimate. It 
is derived for some laws in Appendix A. 

As a conclusion, the volumic and boundary equations of /^[ui, <Ji\ (u, q, 0, cr) = (f , 0, m, 0, 0, 0, 0) 
(see (117)1 for the definition of P 2 ) are satisfied by the limit of the sequence. 

Concerning the initial velocity, u™ converges to u in iJ5+l(0, T; H 1 + r_61 ) with < 6 < r (cf. 
Lemma [3. 5p . So we have proved that u(t = 0) = in H 1+r ~ e and so also in H 1+r . 

Concerning the initial extra stress, cr™ converges to cr in ^^(0, T; H 1+r ). So cr(0) = in 
H i+r 

Up to now we have proved that fT2Tl) has a solution. By lifting the initial conditions, we will 
solve ffl~5l). This is done in the next subsection. 
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3.6 Lift the initial conditions to solve Pi 



To solve the linearized problem about zero deformation (£ = 0), with general rhs, we introduced 
and solved the reduced P 2 problem f[2~Tj) . We want now to solve Pi (or the equation ffTol) ). So let 
(f , a, m, g, k, u , cr ) in a closed ball of radius R in . We also take fi G fT r and mj 6 /T 1+r . 
We start by lifting the initial extra stress. So let (T\ be such that: 



< r 1 + We^l = m 1 

m(o,x) = ff (x) vx 



(41) 



One may easily check that <T\ G L°°(0, T ; H 1+r ) for any T . This and its special regularity 
has been repeatedly used in the estimates because of the nonlinearities. 

Then applying Theorem 4.1 of |4] to the rhs (fi, a, 0, g — <Ti.N, fc, u , 0) provides (ui, q%, 0i, 0) 



in X?p such that 



(42) 



Reu 1>t — (1 — e)Au\ + Vgi 

divui 

-qiN + 2(1 - e)D[u x ] • N - a<9 T (0iN) 

4>i,t ~ <9tUi • N 

ui(t=0) 

The just found ui, ct\ are those denoted the same way as in (l2Tj) . To summarize, these fields 
satisfy 



fi in Q, 

a in J7, 

(/-o-i-N on^ F x(0,T) 

k on Sp x (0, T) 

on S'i? 

onS B x(0,T), 

uq(X) in Q. 



Pi(ui, gi, 0i, oi) = (fi - div a 1 , a, m 1 - 2eD[uJ - We^ (Vui, (Ti),g, k, u , cr ). 

and so (ui + u, q 1 + q, 0i + 0, <j\ + a) G is such that 

P 2 [u 1 ,o- 1 ](u,q,(j),a) +Pi(ui,gi,0i,cr 1 ) = 
Pi(ui + u, qi + g, 0i + 0, o"! + a) = 

(f + fi - div eri,a,m + mi - 2eD[ui] - We^ a (Vui, <Ti),g, k, u , tr ). 

So, by choosing fi = diver! and m = 2eD[ui] + We^ a (Vui, o"i), we provide a solution 
(ui + u, q 1 + q, 0! + 0, o - ! + er) of (fT5|) as announced in Theorem 13.11 Notice that we cannot 
set mi := 2eD[ui] + We5> (Vui, a\) since ui,cri depend on mi. It would be only an implicit 
equation. Then it suffices to let mi be whatsoever and determine m as a function of ui, <j\ as 
done here. 

To complete the proof of Theorem 13.11 we still must prove the inverse of Pi is bounded. 



3.7 Uniqueness and boundedness 

To prove uniqueness of the solutions to ffT5]) and the boundedness of Pf \ we start from the 
equation ([15]) satisfied by (u, q, 0, er) G XJ o for given (f , a, m, g, k, u , cr ) in a closed ball B in 
Yf. First we prove uniform estimates and then boundedness of Pf 1 . 
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3.7.1 Uniform estimates 



To have uniform estimates, we may pass to the limit in (130]) and recover the estimate for the 
velocity 

|(u,g,$,0)Uf o < C, \Vn\ L 2^ T . Hl+r) < C. (43) 

and the extra stress 

Hl+r(t) < S. (44) 

These bounds are true only for initially vanishing fields. But then it suffices to add a lift func- 
tion to have the same bounds for general (u, q, <j), cr) solution of the system ()15p . 

It is not difficult to mimic the proof of the inequality (128]) made in Subsection 13.41 to prove 
even a stronger result. The bound (jHJ) enables to prove that the partial derivative in time of 
cr is also bounded in L 2 (0, T; H l+r ). Indeed, the nonlinearities can be managed thanks to (132]) . 
Similarly to (134]) . it is then easy to get that: 

I CT \H 1 (0,T;H 1 + r )< C (| m |l 2 (0,T;H 1 +'') + I Vu |L2( 0iT;ff l+r) +S) . 

3.7.2 Boundedness 

Let us assume ffT5]) is satisfied by (u, p, (ft, cr) and (u',p', (ft', cr') with two right-hand sides. 



Re<9 t (u - u') - (1 - e)A(u - u') + V(q - q') = f - f + div(<r - cr') 
div(u — u') = a — a' 

cr-cr' + We (d t (cr - cr') - g a (Vu - Vu', cr) - g a (Vu', cr - cr')) = 

= 2eD[u - u'] + m - m' 
(a - a') -N - (q- q')N + 2(1 - e)D[u - u'] ■ N - ad T (((f) - <j/)N) = 

= 9-9' 

(<f) - <j)% - d T (u - u') • N = k - k! 

(0-0')(t = O)=O 

(u-u')(t = 0) = (u -u' )pf) 

(cr-<T')(t = 0) = (<T -Cr' )(X) 

u - u' = 

Theorem 4.1 of |4j gives: 



in Q x (0,T), 
in Q x (0,T), 

in n x (0,T) 

in S F x (0,T), 
in S F x (0,T), 
on Sf, 
in fi, 
in Q, 

on S B V*. 



(45) 



I u-u',q- q',(f>- 4>',0 \ X r T < 

< C I f - f + div(cr - cr'), a - a', 0, g - g' - (cr - cr') ■ ~N , k - k', u - u' 0) \ Y r 

< C \f-f',a-a',0,g- g', k - k', u - u' , \ Y{ + ^ b > 

+ C \(T -cr' | L 2 (0iT; Bi+r )n Bi(0,T;Hi) +C\cr -cr' \ Rr+ i {Sjr><(0jT)) • 

Therefore we must estimate cr—cr' in three norms L 2 (0, T; H 1+r ), H r ^ 2 (0, T; ff 1 ) and if ^ + ^(0, T; iJ 1 / 2 ).! 
We also need estimates of cr — cr' in 77^(0, T; H l+T ) to complete the norm of Xj,. The last 
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norm is the strongest but a priori no bound can be found for initially non-vanishing functions. 
Nevertheless we will establish such bounds because we can lift the initial conditions and then 
use classical results. 

Remark 3.7. By various estimates, we can bound \cr — cr'\ with terms including |cr — cr' \ but 
also |Vu — Vu'|. The latter term (and not the former) needs to have a coefficient sufficiently 
small for T small enough. This will enable these terms to be absorbed by the left-hand sides of 
( Tffi| ) and so lead to the expected bound UiOi) . 

From now on, we assume (u, q, </>, cr) solution of ffT5]) and so initially non- vanishing. If the 
functions vanish initially, the bounds in L 2 (0, T; H 1+r ), H r/2 (0, T; H 1 ) and iJ5+i(0, T; H 1 / 2 ) 
are included in the bound in (0,T; H 1+r ). So the hierarchy of these norms is a priori 
not obvious. Yet we prove below that the if~(0, T; H 1+r ) bound ensures the others. To 
apply Lemma I3.5[ we need to lift the initial values. Then we apply the process that led us 
to (HU1) to the new initially vanishing functions. Moreover, we will even derive an estimate in 
H^O^T; H 1+r ). 
Let r be such that 

r + We d t T = 0, 

T(0) =<7 , 

and in the same manner r'(t, X) = e~vfecr' (X). Instead of (T451U . by renaming cr := cr + r and 
cr' := cr' + r', we have: 

cr - cr' + We (d t (& - cr') - g a (Vu - Vu', <r) - g a (Vu', & - cr')) = 

= 2eB[u - u'] + m - m' + We # a (Vu - Vu' r) + We # a (Vu', r - t') in Q x (0, T), 

(47) 

and mainly (cr — cr')(t = 0) = 0. Moreover, all the uniform bounds stated in Subsubsection 
13.7.11 apply to cr and cr' . Now, we take the H 1+r scalar product of (H71) with cr — cr', then 
simplify by |<r — cr'| 1+r and we have: 

I a — a U +r +We— \ cr — a' |i +r < 2e I Vu — Vu' h +r + I m — m' |i +r + 
at 

+CWe [| Vu - Vu \ 1+r (| cr \ 1+r + | T |i+r)+ I c> - c>' | 1+r | Vu' | 1+r + 
+ I Vu' \i +r \ t-t' |i] . 

This equation writes also, thanks to the fact that |t — r'|i +r (t) < |cr — cr' Q | 1+r , |r|i +r (t) < 
|cr | 1+r and the uniform bound (|44|) : 

(1 - CWe|Vu'|i +r )|<r - a \i +r +We-^-|cr - &' \ 1+r < C(|Vu - Vu'|i +r + |m - m'\ 1+r )+ 

+C|Vu |i +r |cr — cr' |i +r . 

and so 

\cr - &' \ 1+r < |cr - cr^ \ 1+r e~ We So (l-CWe|Vu'| 1+ ,)d S + 

+C/ t e-^^( 1 - CWc ' Vu '' 1+ ^'(|Vu-Vu'| 1+r (t) + 
+ |m - m'|i +r + |Vu'|i +r (t)|cr - cr' |i +r )ds. 
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Since the term in the exponential can be bounded thanks to (143]) : 

J (1 - CWe| Vu'| 1+r )ds' < C J iVu'li+r < Cv^l Vu'| L2(0iT . H i +r) < Cv^, 

we prove in L°°(0, T; # 1+r ): 

| 6- - <x' |i +r (i) < C (| Vu - Vu | L 2(o,t ; /P) + I m - m' |l2(o,t ; hi) + | cr - erg | x ) , (48) 

where C = C(e, B, We, S, V, Re). A similar bound in L 2 (0, T; H 1+r ) is easy to derive. 

The estimate on the time derivative of cr — cr' is similar to the one derived in (140]) (still because 

(cr — <r')(£ = 0) = 0). By taking the H 1+r (Q) norm of the time derivative in (j47|) . we get: 

|<9 t ((X - <X )| 1+r < ^|<T - O-'ll+r + C (|VU - Vu'| 1+r |<x| 1+r + |Vu'| 1+r |<T - &'\ 1+r + 

+ | Vu - Vu'|i +r |r|i +r + |Vu'|i +r |r - r'| 1+r ) + 
+ We' Vu ~ Vu 'li+- + We' m ~ m 'li-H- 

Then we take the L 2 in time norm of this inequality and have: 

\d t (& - o-')| L 2 (0iT . H i+r) < C(|<r o -0"d|i+r + |Vu- Vu'| £ a(o,r;Hi-H-) + |m-m / |£a(o,ri£ri+'-)) J (49) 
and finally we gather all the results in 

\ (T — <t' |.ffl(0,T;.ff 1 +'-)< C (| Vu — Vu' \L 2 (0,T;m+r) + | Ul — Ul' | L 2 (0,T;m+ r ) + \ & Q ~ Ct'q |l+r) • 

Now, since | e \ H s( 0)T )< C for s > 0, we can state the following properties on cr and cr': 
I cr — cr' |L2( 0iT .^i+r) < \t — r'\ L 2^ T . H i+r^ + |cr — a-'\ L 2^ T . H x+r) 

< C | CT — CTq | 1+r +CT £ |<X — & \H 1 (0,T;H 1 + r ) /r n \ 

<- r* I /-r' I _L N 0U / 

Sl> I CTO — CTo U+r + 

+CT e ' (| Vu - Vu' Iz^o^i+r) + | m - m' |z,2 (0iT . H i+r)) • 
The same can be done in if r//2 (0,T; H 1 ): 



\ <T — & \h%(0,T;H 1 ) — l T ~~ T lfl"5(0,T;fl"l) + 1^ ~ ^ l/f^O.T;// 1 ) 

< C | CT — CTq In + CT e |(T — C> Iffl^^if 1 ) 

< C | <T — CTq I ! + 
+ CT €> (| VU - Vu' |£3( 0) T ; jfl+r) + | m 

In order to bound I cr — cr' I r , 1 , i, „, we can do in the same way as above by lifting the 

initial conditions and by replacing the scalar product in H 1 by the one in H 2. We may then use 
Lemma 13.21 and Lemma 13.51 because the initial fields, after lifting the initial conditions, vanish. 
Accurate estimates in if 1 (0,T; if 1 / 2 ) enable: 



(51) 



m ' \L 2 {0,T;H 1 + r )) 
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. cr — cr I r , i , „, N < CT 6 ' ( I Vu — Vu' I , 1 + I m — m' I , 1 ) + , * 

1 l H? + *(0,T;L 2 (S F ))— V l L 2 (0,T;H?) 1 l L 2 {0,T;H? )) (52) 

+ C | CT — CTg I l . 

2 

In order to bound I cr — cr' I i+ r , we lift the non-zero initial conditions and then 

1 ] H~a~ (0,T;L 2 (n))' 

process as above with the L 2 scalar product. Since T is such that Cy/Jo | Vu' \L 2 {o,T ;H 1 + r )< 1 
(cf. f H3|) proved in Subsection I3.7.ip . for T small enough, we have (r < 1/2 < 1) : 



cr — cr' \ 



U ^ {0T .< C I (T - cr' | +CT e ' (| Vu - Vu' |l2(o,t ; l2) + | m - m' |l2 (0iT;L 2)) . 

" (53) 
The inequalities (|50|) . (|53l on cr — cr' can be mixed into: 



cr — cr 



' |jfi+r< C | cr - cr | 1+r +CT e '(| Vu - Vu' |^2+r + | m - m' |l 2 (o,t ; h 1 +<-))- (54) 



Concerning | cr — cr' \ H i^r T H1+r ^ we already have the estimate in L 2 (0, T; H l+r ) (see f l50|) ) and 

L°°(0, T; H l+r ). By lifting the initial conditions, one may, as for (1511) . bound the -ff 1 (0, T; H 1+r ) 
norm and get: 

I cr — cr' I i+r , , , „ < C I cr n — crn 1 1 + r + 

1 l H^I-{0,T;H 1 + r )— 1 " |i+r (55) 

+CT e ' (| Vu - Vu' \ L 2(o,T;Hi+r) + | m - m' \ L 2(o,T-,m+-)) ■ 



We are in position now to use the results above in the intermediate (j4*6j) : 

| u- u',q- q',(f) - 0',O \ X ' T < C | f - f, a - a', 0, g - g', k - k', u - u ,0 \ Y - + , . 
+ CT e ' | u - u' \ K 2+r +C I m - m' \l 2 (o,T;H 1 +^) +C | cr - cr 1 1+r . 

For T sufficiently small (and the constants do not depend on T < T ), the |u — u'|^2+r in the 
rhs can be absorbed by the lhs. So we have: 

| u - u', q - q', (j) - (f)', l^r < C | f - f, a - a', m - m', g - g', k - k', u - u' , cr - cr' \ Y r, 

which may be combined with fl55|) to get the continuity (120]) of P-[ in Theorem I3.ll This 
completes the proof of this Theorem. 

□ 



4 Solving the second auxiliary problem 

We want to solve P<i [ui, cr^u, q, 0, cr) = (f , a, m, g, k, 0, 0) (see ffTTj) for the definition of P2) 
with m, (Ti given and (u, g, $, cr) G Xf = {(u, q, $, cr) e X£/u(0) = , cr(0) = 0}. Indeed, we 
prove the following well-posedness Theorem. 
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Theorem 4.1. Let < r < 1/2, B > be given and (ui,<Xi) in a bounded subset of 



K r+2 x H 2 (0, To; H 1+r ). For any (f, a, m, g } k } 0, 0) m a ball Byf(0, B), there exists a unique 
(u, q, <3>, cr) G Xf solution of 

P 2 [ui,cr 1 ](u, g, $,cr) = (f , a, m, fc, 0, 0). (57) 

Moreover the P 2 [ui, en] -1 operator is bounded and even continuous in Py^(0, B) and its bound- 
edness constant does not depend onT < T nor on u 1; en in a given ballofK^xH^r^^H 1 ^ 

Proof. Incidentally, while solving the first auxiliary problem, we found a solution of (J2TJ) which 
happens to have the same lhs with some vanishing rhs. So we are going to lift the rhs of 
conditions f|T7|) 9. ffT7|) ^. f fTT|k so as to solve P 2 [u 1; en](u, g, 0, cr) = (f, a , m , 9, k, 0, 0) which is 
the same as floTl) . 

It is proved in [4J that if (0, a, 0, g, k, 0, 0) G , there exists (u 2 , g 2 , 02, 0) G X£ with continuous 
dependence such that: 

' Re u 2 , t - (1 - e)Au 2 + Vg 2 =0 

divu 2 = a 

-g 2 N + 2(l-e)D[u 2 ] -N-a<9 T (0 2 N) =g 

2 , i -(9 T (u 2 ))-N = fc 

[ u 2 {t = 0) =0. 

Instead of solving (l2~li) . we could have solved the same system with ui replaced by ui + u 2 and 
with an other rhs: 



infix (0,T), 
in VI x (0,T), 
on S F x (0,T) 
on S F x (0, T) 



P 2 [ui + u 2 ,cn](u,g,0,cr) = (f , 0, m + 2eD[u 2 ] + We 9a (Vu 2 , en), 0, 0, 0, 0). 
Such an equation can be rewritten 

P 2 [u 1 ,cr 1 }(u,q,<P, ( T) = (f,0,m + 2 £ D[u 2 ] + We{g a {Vu 2 , cr) + </ a (Vu 2 , <n)) , 0,0, 0,0). 
Let us write explicitely the Stokes part of this system: 
Re d t u — (1 — e)Au + Vg — diver 



div u 

cr • N - gN + 2(1 - e)D[u] • N - a<9 T (</>N) 
t - (&ru) • N 

<f>{t = 0) 
u(t = 0) 
u 

that can be rewritten thanks to the definition of (u 2 , g 2 , 2 ): 

Re d t (u 2 + u) - (1 - £)A(u 2 + u) + V(g 2 + g) - div(cr) 
div(u 2 + u) 

cr • N - (g 2 + g)N + 2(1 - e)B[u 2 + u] • N — ad T ((<p 2 + ■ 
(0 2 + 0)t-(<9 T (u 2 + u))-N 

(02 + 0)(* = O) 

(u 2 + u)(t = 0) 
u 2 + u 



f infix(0,T) 

infix(0,T) 

onS F x (0, T) 

onS F x (0, T) 

in tt, 

in fi, 

onS B x (0, T) 



= f infix(0,T) 

= a infix(0,T) 

)N) = p onSpx(0,T), 

= k on5 F x(0,T), 

= in fi, 

= in fi, 

= onS B x(0,T). 



(5? 
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Now we must write the viscoelastic part of the system 



cr + We (d t tr - ff a (Vu, cr) - ^(Vuj, a) - g a {Vu, <xi)) - 2eB[u] = 

= m + 2eB[u 2 ] + We (g a (Vu 2 , cr) + g a {Vu 2 , (t 1 )) infix (0, T), 

that can be rewritten 

o-+We (d t <r - g a (Vu 2 + Vu, cr) - g a (Vu u cr) - g a (Vu 2 + Vu, cr 1 ))-2eD[u 2 +u] = m. (59) 
Equations (158]) and (|59|) can be rewritten 

P 2 [ui,cr 1 ](u 2 + u, q 2 + g, 02 + 0, cr) = (f, a, m,#, A;, 0,0), 

which is exactly (|57|) in Theorem 14.11 for (u 2 + u, q 2 + g, 02 + 0, cr) £ -^f- 

Uniqueness of the solution and boundedness of this operator P 2 _1 [u 1 , cr^ may be proved in the 
same way as for the operator P^ 1 and even its continuity (it is nonlinear). The proof relies on 
the boundedness of ui, <T\ in the appropriate spaces. 

This completes the proof of Theorem 14.11 

□ 



5 Estimates of the error terms 

In fIBUl) we expand P in the neighborhood of £ = 0. If we keep only zeroth order terms of £ 
for P(0, 0, 0, 0, 0) and Pi, we put the "error" terms in E. Thanks to the fact that is initially 
vanishing (which is not the case for the other fields like the velocity), we linearize in 0, put 
apart zeroth order terms in P(0, 0, 0, 0, 0), put the linear (in 0) term in Pi(u, q, 0, cr) and higher 
order terms in E. The operator Pi is nonlinear because of our nonlinear constitutive law. So, 
one may write the full system of equations to be solved: 

Order 0(£)+0(fli 2 ) 

P(£,u,g,0,<r) = P(0, 0*0, 0, 0) +Pi(u, q, 0, cr) + P(£, u,g, 0, c^ (60) 
= (0,0,0,0,0,u ,cr ), 
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where the E terms are fully written below with usual summation convention and = 

1 - e)dj{iijUi,i) + ikidkq - <Jij,h£,kj 



E}{£, u, q,<f),cr) 



Zkj u j,k 



M(£,u,g,0 )C r) = -We 



2 ~ (&i U k,l°~kj + OikUkjilj) + 



(a+ 1 



H 2 ( a ikilkUjJ + Ui^lk&kj) 



^ 4 (£,u,g,0,cr) 

where Qi 

and Q2 
E 5 (£,u,q,(/), a) 

£ 6 (£,u,g,0,cr) 



- e(u^kikj + u jt ki 



ki 1 



-q{Mi - iVj) + (1 -_e) [^fcj«<,fc + 6wUj,*] 
+(1 - e) [£ kj u ijk + £ ki u jik ] {Nj - Nj)+ 
+g((X 1 )(Af i - Ni) + - ad T Q { + a^Mj - Nj) 



(61) 



'1 + h 



2\+- 



1 + 



-1 + 



cv 



1 + h 



l + ti 2 

-('^(i + h'^ + tt + OQi 
-d T (u) ■ N [(N 2 + d Tr]l y 2 - N 2 2 ] - 

(N 2 + <9 T ni)~ 2 (d T u 2 d T ri 1 - d T uid T 7] 2 ) 
O = £ 7 (£,u,g,0,cr). 



Lastly, we will lift the initial conditions with (u 1; q x , 1; o^) in . We will also use the space 
of initially vanishing fields X? = {(u, q, 0, cr) G X^/u(0) = 0, <r(0) = 0}. We will denote 

(0, -R) the ball of center and of radius R in 
Actually we need to bound the "error" terms. This is done in the following theorem. 

Theorem 5.1. Let < r < 1/2, and (ui, gi, (fix, cr x ) G Bjqr (0,-R). There exists e' > and 
< Tq < T depending on (iii, gi, 0i, <Ti) and i?, sncn tnat if < T < Tq and (u,g,0, cr) G 
(0, -R), i/ien i?(ui + u, gi + g, 0i + 0, + cr) is in i/ie space YJ(f2), and t/ie following holds: 



I £ 4 ( Ul + u, 9l + g, 0! + 0, a x + cr) | (y .). < CT £ '. 
//in addition (u', g', 0', <x') G Bx*{0,R), the operator E is contracting: 

I £(ui + u, gi + g, 0i + 0, o"! + cr) - £?(ui + u', gi + g', 0i + 0', cr x + cr') | K r 



< CT e ' I u - u', g - g', . 



(62) 



(63) 



w;i£/j a constant C that depends on e, a, We, r, R, (u x , q x , 0i, <Ti) ; 6nt not on T provided T < Tq. 

We provide the proof for the first and third components more specific to viscoelastic fluids. In 
the fourth component {E 4 ), we study only the term added by viscoelasticity of the fluid. The 
other components (E 2 , E 4 , E 5 ) are studied in [I] or in [B]. 



5.1 Various lemmas 

We want to give a meaning to the fact that £ is small in small time (used in Theorem 15. ip . So 
we need some lemmas stated hereafter. 
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Lemma 5.2. Let < T < T , < s < < e' < s and X be a Hilbert space. The linear 

mapping v ^ V(t) = / v(s) ds, is a bounded operator from H s (0,T; X) to H s+1 ~ £ '(0, T; X) 

Jo 

and 

I V | a+w < CT 6 ' | v | s , 

with C independent of < T < T . 

Proof. The proof is through double interpolation. For s = 0, Cauchy-Schwarz inequality gives; 

| V \ X if) < t* | V \h°{0,T;X) =H ^ |o< I U |o, 

where |.| = \-\h°(o,t-,x)- So | V |i< (1 + CT 2 ) | u | < C | v | where |.|i = \.\ h ^(o,t-x)- Then, by 
a classical interpolation inequality (1 — e' > 0), the result is obvious for s = 0: 

I V | w < C | V U^'l V |£< CT e ' | v | . 
In the same way if v G i? 1 (0, T) 

I V li = l v \l + | V \l< C | u |g if T<T . 

Still if t> G 1 : | V ]2< C | v |i. Moreover if v(0) — (this assumption is discussed below), 
since i> is the integral of d t v, then | t> | < CT \ v' | and |V|i < CT|i>|i. The same interpolation 
inequality provides: 

I V \ 2 -e'< C | V V \(< CT 6 ' I v \(\ v |}- e '< CT e ' | v |i . 

By re-interpolating between the two inequalities, one completes the proof. 
Notice that the assumption v(0) = disappears since we restrict the regularity to s < 1/2. 
Such an assumption is meaningless for non-regular functions. So this assumption, also done by 
J.T. Beale (in his Lemma 2.4), does not limit the proof. 

□ 

Let us remind that we define -ff _1 (fi) as the dual space of °H 1 (Q) = {pG if 1 (f2), p = on Sf}- 
This space is needed for the incompressibility condition. We may then state our next lemma. 

Lemma 5.3. Let Q be an open subset o/M 2 . 

(i) Let r > l,r > s > ; v G H r (Q) and w G H S (Q), then vw G H s (£l) and 

| VW \ s < C I V \ r \ W \ s 

(ii) Let v G H r (Q) for r > 1, and w G oH^^Q), then vw G oH' 1 ^) and : 

I vw |_i< C | f | r | w |_i 

fiiij Lei v,w G if 1 (f2), £/ien utti G L 2 (f2) and 

I fiu |o< C | f |i| w |i 
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(iv) Let v G H 1 ^) and w G L 2 (Q), then vw G qH l (Q) and 

I vw |_i< C | v |i| w |o 

This lemma is proved in dimension 3 by J.T. Beale (Lemma 2.5 p. 366 of [6]). The proof in 
2-D is very similar. Following J.T. Beale, we state below that the product of two functions in 
appropriate spaces is a continuous map. 

Lemma 5.4. Let X, Y, Z denote three Hilbert spaces and M : X x Y — > Z , a bounded and 
bilinear map ("multiplication"). 

(i) Suppose u G H s (0,T;X) and v G H 8 (0,T;Y) where s > 1/2 then uv = M(u,v) G 
H s (0, T; Z) and \ uv \ s < C \ u \ s \ v \ s . 

(ii) If s < 2 and u,v satisfy in addition to (i) the conditions <9 t fc n(0) = = d^v(0) for 
< k < s — 1/2 and s — 1/2 is not an integer. Then the constant C of (i) does not depend 
on T < T . 

This lemma is proved p. 366 of [Ej (Lemma 2.6) in two lines. The proof is based on an extension 
to t G M and estimates of the transforms of (i). Lemma 13751 enables to prove (ii). 

We will also need the following lemma that ensures there exists a bound independent of T < T 
for fields whose initial value is not zero. 

Lemma 5.5. Let (ui, q%, 4>i, o*x) G XL (r < 1/2) and (u, q, 0, cr) in a ball of radius R in X^ 
with u(0) = 0, cr(0) = (so in X^,). Then the following holds with constants C that depend on 
R, Ui, qi, 0i, (Ti, r, not on T < T : 

| fc c?j(ui + u) |xT(o,r ; n)< C (64) 
| 9 fc (ui + u) |ifr + i (0)T; n)< C (65) 
I V(gi + g) |jfr( 0) T ; n)< C (66) 

| (<Ti + 0")y,fe |jfr(o,T;n)< C (67) 

The proof relies on two arguments. On the one hand (ui, qi, <p±, 0"i) does not need to be 
extended until T . On the other hand u, q, <p, cr are either initially vanishing (for u and cr) and 
so can be extended on (0, To) with a bound independent of T < T , or not very regular (p and 
0) and so the extension on (0, T ) does not generate any cost nor need any extra assumption. 
At that level, we see that improving to r > 1/2 our result would modify the whole proof. 

Proof. Concerning f[64"j) . Lemma 2.1 of [6] (which deals with the boundedness of the derivative 
and trace operators) and our Lemma [3721 (which deals with the extension operator from (0,T) 
to (0, +oo)) enable to state: 

| d k dj(ui + u) \ K r(o,T;n) < | dkdjUx |^( ,t o; q) + I 9 k djU \ K r(o,T-,n) 

< C | Ui |_ftT''+ 2 (0,To;^) +C | U |£,2( ,T;H3+r) f] H r / 2 (0,T;H 2 ) j 

< C | Ui \K r + 2 (0,To;U) +C | U \K r + 2 (0,T ;n), 

< c. 
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One may prove in the same way fl65l) . 

Concerning the pressure, since r < 1/2, we do not need to force unphysical vanishing initial 
value to have the result thanks to Lemma 2.1 of [6]. 
The extra stress bound can be written: 



| (<Ti +Cr)j |jfr<| (Txj \K r (0,T ) + | fj \Kr(0,T)< C ' (T ) \ (7 \ \K r + 1 {Q,T ) +C \ (J |kH-1(o,T )< C. 

□ 

The following lemma gives a sense to the claim that "£ is small" . 

Lemma 5.6. Let r be such that < r < 1/2, and let us denote A the algebra 

A = H 1+ i(0,T;H l+ $(Q)). (68) 

For any Ul G K 2+r (0, T ), u, u' G K 2+r (0, T), u(t = 0) = = u'(t = 0), and £ defined by 
£(u) = (I + drf)^ 1 — I = (I + J Q Vu) _1 — I, there exists e' > and constants C such that if 
T < T ; 

l^ + ^l^^^^l^ + ^U < CT" (69) 
| £( Ul + u) - ^(U! + u') | H i+r (0)T; ^ 1+r/2 <| C(u x + u) - ^( Ul + u') \ A < CT* \ u - u> 

with C, dependent on u 1; r, T but not on T provided T <T . 

The following Remark must be taken into account before we prove this lemma. 

Remark 5.7. The continuity bound of the product in the algebra A depends onT < Tq if there 



is no more condition. If we add the condition that the fields are vanishing initially, Lemma 5^4 



(proved with an extension operator whose properties are given in Lemma \3. fy) may ensure that 
the bound does not depend on T < T , but what happens if this assumption is not satisfied ? 
To answer this question, let us write the continuity of the product of two functions that do not 
depend on time and so do not vanish at t = 0: 



1 x 1 ltf 1+ §(o,T)^ G ( A ) I 1 U 1+ S(o,t)I 1 \n 1+ i(o,T)^ 1 ^ C ( A ) I 1 ltf 1+ £ 



(0,T) • 



So the constant not only depends on T, but even tends to +oo when T — > 0. 
One might wonder whether the Lemma lSTR needed in the proof of Lemma \5Jh uses a too specific 
extension which could be improved. Of course, the extension used will not work for instance on 
constant functions which will not remain in any H s (0, oo). But the above inequality proves that 
no other extension operator could suit. This explains why in Remark \3.!% we claimed we were 
forced to lift initial conditions so as to have new initially vanishing fields. 

Proof. We denote £ := £(ui+u), £' := £(ui+u') and rj, rj' the associated fields respectively. We 
repeatedly use the fact that A is an algebra and that u is in K 2+r and so also in H~ (0, T; H 1 ). 
So £ G H\0, T; H 1+r ) f| H^(0, T; L 2 ). 
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Let us denote S the extension to (0, +00) operator and 11 the restriction to (0, T) operator. 
From Lemma 13721 we know that S is bounded independently of T < To for convenient fields. So: 

I drjdrj' | H i+g (0T .^ 1+r/2) =| US drillS drf \ (0)T;iJ - 1+r/2) 

— ^ I 7t£& r ) l# 1 +s(o,T;H 1 +'V 2 )l H-Sdrf |^i + g (0T . Hl+r/2 ^ 

< C I Sdrj l^i+i^.jyi+r/a)! £dr/' | Hl+ g (R . Hl+r/2) 

< C I drj U| dry' U, 

where C does not depend on T < T . As a consequence there exists cti such that if | dr] |^< oti, 
then, £ = (1 + drj)- 1 - I is in A and it can be expanded in series with: 

I £ U< C I dr, \ A . (71) 

Since dr?(X, t) = / V(ui + u), u x e K r+2 (0, T ) and the fact that u is only in K r+2 (0 } T) but 
Jo 

with vanishing initial condition, Lemmas 15.51 and 15.21 enable to state the announced result f|69jl : 



I I U< C I drj | H i+r/8 (0iT;H i+r/2)< CT e | V(ui + u) | H &(o,T; J ff 1 +'72)- • 
Moreover, since *4 is an algebra, if | dr? |< «i/2 and the same for dr/'; 

£ — £'=(! + £)d(r/ - r/)(I + £') = d(r, - r/) + £d(r, - r/) + d(r, - r,')g + £d(r7 - r,%' . 

Here, all the functions are initially vanishing. So, in the same way as for flTTj) : 

I \a< C I d(r, - r,') U • (72) 



One may then prove (1701) in the same way as fl69l) . □ 

The previous lemma will not be sufficient in some estimates. So we state the following lemma 
which will turn to be useful. 

Lemma 5.8. If the velocity u is in L 2 (0, T; H 2+r ), then the field 

£ = (Id + dr/) 1 - Id = (Id + f Vu) 1 - Id 

Jo 

is in H\0,T;H 1+r ) with 

I £ \H 1 (0,T;H 1 + r )^ C. 

Moreover, for < r < 1, we define the algebra 

A' = H^{0,T;H 1+r {n)), (73) 

for which 

\$U><CT*'. 
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The proof is similar to the one of Lemma T5.6I for the first estimate. Then one must use Lemma 
13.41 to bound \ £\a' with a CT e and | £ \m{Q,T-H^+ r )- The rest is left to the reader. 
In the proof of estimates in L 2 (0, T; H s ) below, we will use Lf 2 estimates and then H ^ (0, T) 
L°°(0, T). The constants in Sobolev's inequality must be independent of T < T . In the general 
case, it is wrong but for the subspace of initially vanishing functions, it is true as states the 
following lemma. 

Lemma 5.9. Let < r < 1. If v E H~(Q,T) and is initially vanishing then 

I v |l» (0jT) < C I v \ H i+r^ T) , (74) 

and the constant C does not depend on T < T . 

Proof. The proof of L°°(R) is classical. Then, thanks to the properties of the 

extension operator E in H l (Lemma 13.21 (ii)), 

I V |l°c(0,T)<| £(v) |x°°(K)<| \m(M)< C I V |/fi(o,T), 

with C independent of T < T . 

The same can be proved for non-integer Sobolev spaces and completes the proof. □ 

For the estimates in H%(0, T; L 2 ), we will need a more precise result than Lemma 15^41 to estimate 
the product of two functions. 

Lemma 5.10. Let X,Y,Z denote three Hilbert spaces and M : X x Y — > Z , a bounded and 
bilinear map ("multiplication" ). Let 1/2 < s < 3/2 and < s' < s. 

(i) Suppose u e H s {0,T;X) and v G H s '{0,T;Y) then uv := M(u,v) G H s '{0,T;Z) and 
| UV | s '< C I u | s I V \ s /. 

(ii) Let it, v satisfy the conditions of (i) and in addition u(t = 0) = and s' < 1/2. Then 
there exists a constant C in (i) that does not depend on T < T . 

Proof. Let us consider the functional v i— ► uv. For u G H s (0, T; X), this function can be defined 
in H°(0,T; Y) and in H s (0,T; Y). It satisfies: 

I UV \hO(0,T;Z)< C I U \h°(0,T;X) \ V \h°(0,T;Y), 
I UV \h s (0,T;Z)< C I u \h s {0,T;X) \ V \h s (0,T;Y) ■ 

Then a simple interpolation provides the result (z) for s' between and s. 
Then the extension operator enables to exhibit C independent of T < T : 

I uv I/p'(o,t)<I S{u)S{v) \ H s> m < C | £(«) ^(v) Ih-'cR)- C I u l^ a (o,r)l u \h*'(o,t) ■ 

□ 
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5.2 Estimates on E 

The terms to be estimated are in flfJTj) . 

For the Navier-Stokes equations, G. Allain [I] sends back to [3] where she indicates the tools 
to get accurate estimates and refers to J.T. Beale's article [5] for details. 

We denote £;y(ui + u) = £jy = ((Id + drj)^ 1 — ld) k j because it depends on the velocity ui + u 
through r). For the sake of completeness we consider in detail the first term of E l : £jy(ui + 
u)d fc [Cy(ui + u) x (m + u) M ]. We prove it is bounded in K r (0,T) by CT e> (see Qe5)) and 
contracting (cf fl63|) ). 

l + r 

The main difficulty here is that the constants for the embedding (H~(Q,T) L°°(0,T)) or 
the property of algebra (| |h s < C | / | i±r \ 9 \u s for < s < 1/2) tend to infinity when T 
tends to in the general case. Yet, since the fields involved are initially vanishing, we can use 
our more refined lemmas proved above. 

Moreover, for the reader not familiar with the K r spaces, we prefer to split the estimates in 

L 2 (0,T;H r ) and H§(0, T; L 2 (f2)). 

Below, all the constants C are independent of T < To. 

In L 2 (0,T; H r ). 

Here, we can use the algebra L°°(0,T; H 1+r ). Then, denoting £ = £(ui + u), 

I £kjdk[€lj X ( u l + u )i,l\ \L 2 (0,T;H r ) — C \ £ kj \L^(0,T;H 1 + r ) I ^fcjCzj X ( u l + u )i,l] U a (0,T;H r ) 

< C | i kj | H l+r (oT; ^ 1+r) | X (Ui + u) i(J \mo,T;m+r), 

thanks to Lemma 15.91 Then, one may use Lemma 15.81 and write £ = Id + £ to pursue the 
bound: 

< CT e ' (| (Ui + u)i,j) \L 2 (0,T;Hl+r) + \ Clj x ("l + u) M ) |l2(0,T ; H 1 +'')) • 

Then one has, with arguments similar to above and with £y(t = 0) = 0: 

< CT e (| Ui + U |L2( 0)T;H 2+r) +CT e I Ui + U \l 2 (0,T;H 2 + t )) 

< CT e ', 

where C does depend on R but not on T < T . 
In gi(0,r; L 2 ). 

We can no more take off the £ term by a simple L°° bound. So we need algebra properties (in 
time) and so Lemma 15.101 which is a more refined version of Lemma 15.41 Indeed the latter is 
useless in H^(0, T; L 2 ) since r/2 < 1/2 and so H r l 2 is not an algebra. 
Because of Lemma 15.101 and A = H 1+ s (0, T; H 1+ z) , 

I Ckjd k [lij x (ui + u) M ] | H 5 (0jT . L2) < C | i kj \a\ dh[Iij x (ui + u) M ] 1^/3(0,^) 

The last inequality uses Lemma 15.61 Then a simple decomposition of £ = Id + £ enables to 
have a (ui + u)u term not initially vanishing but linear, and a £(ui + u)(ui + u)^ term which 
is nonlinear but initially vanishing: 
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< CT e (| Ui + U \H r / 2 (0,T:H 2 ) +CT e \ 111 + U \h t / 2 {0,T;H 2 )) ■ 

The difference with the estimate in L 2 (0, T; H r ) is that the equivalent term was | ui+u |l 2 (o,t ; h 2 +'-)-| 
Such a term could be bounded by | ui + u \k 2 +>-(o,t) with a constant equal to 1. But here we 
need to use Lemma [3751 and especially its (ii) because we need constants independent of T < T . 
So we split U! + u, use the fact that u x e K 2+r (0, T ) (costless) and that u(t = 0) = 0, so that 
u G K 2+r (0, T) can be extended without any loss to K 2+r (0, T ). In summary, we have proved 
the following bound: 

| £iti(ui + u)<9 fe [^-(ui + u) x (m + u) M ] l^s^.^) < CT € '. 
The other terms are no more difficult. So (1^2]) is established for E 1 . 

Remark 5.11. The argument that t](ui + u) = J *(ui + u) is initially vanishing and so that 
the embedding of in L°° has a constant independent ofT<T is not in p. 380 of [6]. We 

only added this and Lemmas 15. ,91 and \5.10\ to complete the proof. 

One could prove in the same way that E 1 is lipschitz and so satisfies fl63|) . 
5.3 Estimates on 

We consider now the error terms of the constitutive equation in K r+1 ((0,T) x Q) : 



| J5?.(£(ui + u), ui + u, qi + q, 0i + 0, cr 1 + cr) \ K r+^< 
C | £/i(ui + u) x (ui + u)^(<ri + cr) kj \ K r+i +C | (ui + u) itk £ kj (ui + u) \ K r+i . 

In order to treat the first term in K 1+r , we split K 1+r in L 2 (0, T; H 1+r ) and #^(0, T; L 2 (ty). 
We apply Lemma [5.51 Lemma [5.81 and Lemma [5.91 to the first term in L 2 (0,T; H 1+r ): 

I &t(ui + u)(ui + u) A . jZ (cr 1 + er) kj |l 2 (o,t ; h 1 +'-) 

<| 6l( U l + U ) U'l ( U l + U )k,l |L 2 ( 0iT;// l + r)| ((T 1 + CT) kj |Loo( 0jT;H l + r) 

< CT £ ', 

where A' is defined in ( 1731) . One may also bound the first term in if 1 ? (0, T; L 2 ). Indeed, 
Lemma [5.61 and Lemma [5.1UI enable to write: 

| £ M (ui + u)(ui + u) fej i(<ri + a-) kj \ H ^ (0iT . L2) 

< C | f K (ui + u) \ A , x 

X(| ( Ul ) M ^l+r (0)To;i2) + I Uk,l \ H l^ {m2) ) \<Ti + <T \ H ^z (0tT , m+r) 

KCT'Wu^qu^ar \x h +CR) 2 

< CT e ', 

with C depending on We, a, ui, qi, <p\, cr 1; T , but not on T < T . So we proved: 
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I £ K (u! + u)(ui + u) M (<Ti + cr) kj \ K r+i< CT e '. 
In order to treat the second term, the algebra properties of A' = {§,T\ H l+r ) give: 

| (m + u)i^^(ui + u) |ifr+i< CT 6 ' | ui + u \ K ^< CT 6 ' , 
thanks to Lemma [5. 51 

In order to prove the contracting property, we compute the difference: 

| £j(£(ui + u), ui + u, q 1 + q, 0i + 0, <ti + a)- 
£j(£(ui + u'), Ux + u', qi + q', 0x + 0', a 1 + a') \ K r+i< 
AC | (f K (ui + u) - £ K (m + u'))(ui + u) fc)J (o-i + cr) kj \ K r+i + 
+4C | f K (ui + u)(u - u') fc){ (o-i + (r) kj \ K r+i + 
+4C | f M (in + u)(ui + u') kj i((T - a') kj \ K r+i + 
+2C | (u - u') i)fc £ fej (ui + u) \ K r+l + 
+2C | (ui + u')j )fc (^j(ui + u) - ^ fc j(ui + u')) \ K r+i . 

Every term can be managed in the same way as before by using (ITUl) and not fl69|) of Lemma 
15.61 Lemma [5 .S\ Lemma [5.91 and Lemma [5. 101 



5.4 Estimates on E A 

This operator comes almost only from the Navier-Stokes part. This part has been estimated by 
G. Allain [3] for a Newtonian fluid. So we only have to estimate the terms (cri + <r)ij(J\fj — Nj) 
in K r+ ^{Sp x (0, T)). We will estimate separately cr 1 + cr and Aj — iVy. 
Concerning <T\ + cr, as for Lemma [5.51 we have: 

| <x 1 + C r \ kT+ x {S F x (0,T)) < | m I^ (5fx(0)T)) + I <r \ K r + ^ SFX{0iT)) 



_ i 



- I CTl lx^(.S F x(0,To)) +C< I ^ l K r +?(S F x(0,T )) 

< | Ui, q U 01, O"! | X r +C | u, q, 0, CT \ X r 

< I 111,41,01,0"! \ X r +CR, 



thanks to Lemma l3~31 (n). So cri + cr is bounded in K r+ ^{Sp x (0, T)) by a constant independent 
of T < To. Concerning Afj — Nj, we use the following formula: 

Af-N= / f (-^T((ui) 2 + (u) 2 )( S ),9 T ((ui)i + (u)i)( S )) ds. 
Jo 

The estimate of J\fj — Nj is done in the algebra H l+ i (0, T; if 5+5 (SV)) so as to conclude. Since 
<9t(ui +u) G L 2 (0,T; ifa+i (Sp)), Lemmas 15.21 and l5~5l apply and also classical theorems found 
in EL8]: 
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L 2 (0,T;H? + %(S F )) 



| J (-a T ((ui) 2 + (u) 2 )(s),a T ((ui)i + (u)i)(s))ds \ H i+i {0tT . H $+t {SF)) 

< CT*' I (-9t(( Ui ) 2 + (u) 2 ), MMi + (u)0) 

< CT e ' | V(ui + u) 1^(0^1+5) 

< CT £ '. 

Thanks to the fact that l + §>2>2[ + 5> ^ ne term satisfies: 

,%3K 3 31 l K r+ ?{S F x(0,T))— 1 1 lK 1 ( b F> 1 3 3 l H 1+ S(0,T;H? + ?(S F )) 

< CTt'l^ + trl^^KCT*'. 
The contracting property of this operator is proved in the same way. 



This completes the proof of Theorem 15.11 

6 Fixed point 

We want to solve the Lagrangian nonlinear system (I14p associated to the operator P. Let us 
remind the reader with the expansion done in the first Section: 

P{£,u,q,<t>,<r) = P{0,0,0,0,0) + P 1 {u,q,<t>, < T)+E{&u,q,<t>,(T) 

= (0,0,0,0,0,u ,<t ) 1 ) 

where P(0, 0, 0, 0, 0) = (0, 0, 0, g ((Xi)Ni - ad T (T), 0, 0, 0) contains all the zeroth order terms 
(gravity and initial surface tension). 

First we lift the initial conditions and the zeroth order terms. To that end, we use Theorem 
13.11 that states that P\ is invertible from Xj, to Yf: Pi(u, q, 0, cr) = (f , a, m, g, k, u , er ) with 
continuous dependence on the initial conditions. So let (u 1; qi, <px, <Xi) G X^ q be such that: 

Pi(ui, qx, 0i, * x ) = (0, 0, 0, 0, 0, u , <r ) - P(0, 0, 0, 0, 0), 

which is allowed thanks to the fact that the right-hand side is in K£ o (£ G H^ +r implies that 
<9 T (T) G P r+ 5). 

If we perform a change of variable for the unknown fields (ui + u replaces u and so on), we are 
led to find (u, q, 0, <x) G X£ — {(u, q, 0, cr) G X^/u(0) = 0, <r(0) = 0} such that: 

Pi(ui + u, qx + q, 0i + 0, o"! + cr) + P(£(ili + u), u : + u, q x + g, X + 0, + cr) = 

P^U!,?!^!,^) 

P 2 [ui,cr 1 ](u, q, 0, cr) = -P(£(ui + u),Ui + u, qx + q,4>x + (fi,cr 1 + cr), 

where P 2 is the second auxiliary problem introduced above in ffTTl) . Since P 2 [ui, crx] is invertible 
when the rhs is in thanks to Theorem 14.11 we want to solve 
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(u,q,(f),cr) = P 2 [ui,o-i](-E(^(ui + u),Ui+u, gi + g, + crx + o-)) 
:= F(u,g,0,cr). 

At that level, G. Allain proves in [3] for a Newtonian fluid that her F lets an invariant ball. 
But since we use a contraction mapping, we do not need this. 

We need now to prove that F is contracting. Thanks to Lemma [5.61 and especially (I63p : 

| F(u, q, 0, a) - F(u', g', 0', a') \ x . < CT e ' | u - u', q - q', - 0', <r - tr' U* , 

where C depends on the parameters, and also on the functions that lift the intial conditions, 
but only through a bound of their norm. As before, it does not depend on T < T , so we have 
proved that for T sufficiently small, F is contracting. We may apply the contraction mapping 
principle which provides a solution to the Lagrangian system of equations fll4p . Owing to the 
regularity of the Lagrangian velocity u, one has solved also the Eulerian equations by a 
simple change of variables. 

Thanks to the contraction property, uniqueness is obvious. 

Last the solution depends on the initial conditions directly in a continuous way and indirectly 
through the lift functions also continuously. 
The proof of our main Theorem 12.11 is complete. 

A Appendix 

The proof of our main theorem is written for the specific constitutive law of Oldroyd. In 

order to convince the reader that it can include most reasonnable constitutive laws, we study 

successively various models well-known in the literature. We could even treat a new model of 

an elastoviscoplastic fluid [21], but it would require deeper modifications of the proof. 

The most crucial step in the proof is the first one: when we prove uniform (in n and T ) estimate 

on a n . Then, using this estimate, it is easy to derive the same inequalities as ours until the 

conclusion. 

A.l The Giesekus constitutive law 

In [11] , Giesekus provides a new constitutive law: 

F) [vl 

a + We ^T <T + C GiesekusV 2 = 2eD[u], 

where We, e, cciesekus are positive parameters. We propose to discretize a reduced version of it 
in which we keep only the new terms: 

* n+1 + Wed t <r n+1 + c Giesekus a n <r n+1 = 2sB[n n ]. 
If we take the H 1+r scalar product of this equation with cr n+1 , one has 

|_n+l|2 | We d /i n+l|2 \ <r (O c \T7-,, n \ _l_ r. It"! \\* n + x \ 

I* 7 1 1+r ~2~&t (/£|VU |i+ r + CGiesekus\<T \l+r)\<T 1 1+r • 
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The same computations as in our proof leads to 
I " n+1 \l + r (t) < 



l+r 

/ (1/2 — CGiesefcusWe | CT n | 1+r ) , 



C (e Js We | Vu n |? +r )ds. 

And as above, we may bound the term in the exponential: 



O f (I/ 2 CdesekusWe | <T n \ i+r ) 2c G iesekus f* l-.nl 

_2 y s ^ dt - We - J s |cr |l+rd * 

- We 

So if we require that T be such that 2cGi eS ekusST /We < 1, we have 

K +1 |i +r (t) < c(|Vu n u 2(0iT; ^)) < cv. 

It suffice then to require CV < S to get the uniform in n estimate on |cr n+1 
A. 2 The Phan-Thien Tanner constitutive law 



l+r- 



In [21], Phan-Thien and Tanner derive new models of viscoelasticity from a molecular argument 
using a network which is allowed to be non-affine: 

Y epTT (trtr)tT + We^U = 2eD[u], 

where Y epTT (x) = exp (epTyWex) in the exponential model and Y(x) = 1 + EpTrWex in the 
linear model. 

We propose to discretize the exponential law in which already treated terms are removed: 



n+l 



+ Wed t (T n+1 = 2eB[u n ] + (1 - e £pTTtran )(T n+1 . 



We take the H l+r scalar product of this equation with cr n+1 and make the same computation 
as above to have: 



<T n+1 | l+r (*) < 



* (\ — M _ gEpTTtro-* 1 



l+r J jj.ii 



t V ' ^dt 



C I (e J if We (| Vu" \ 1+r ))dt'. 

Jo 

One is led to estimate the integral: 

ft (i I i pSpTTtrcr"- I \ 

-/ 1 1 l ^-dt" < C f!, I \ dt » 

J t , We Jt 

< CX*|ES4rr(tr^)*| 1+r dt" 

< C{e £pTTS -1)T . 

So a condition on the PTT parameter Sptt enables to ensure the uniform in n (and in To) 
estimate on |cr n+1 | 1+r . 
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